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Abstract. The paper studies a natural n-dimensional generalization of the 
classical nonholonomic Chaplygin sphere problem. We prove that for a spe- 
cific choice of the inertia operator, the restriction of the generalized problem 
onto zero value of the SO(n-l)-momentum mapping becomes an integrable 
Hamiltonian system after an appropriate time reparametrization. 



Nonholomic systems are not Hamiltonian. Apparently, Chaplygin was one of the 
first who considered a time reparametrization in order to transform nonholonomic 
systems to the Hamiltonian form 9 . Also, after [H], one of the most famous solvable 
problems in nonholonomic mechanics, describing the rolling without slipping of a 
balanced ball over a horizontal surface, is referred as the Chaplygin sphere, see 
[11 [19]. It is interesting that the Hamiltonization of the system by the use of a 
time reparametrization was done just recently by Borisov and Mamaev [U [S] (for 
a geometrical setting within a framework of almost Poisson brackets, see [21)). 

Fedorov and Kozlov constructed natural n-dimensional model of the Chaplygin- 
sphere problem and found an invariant measure [14j . Various aspects of the problem 
are studied in [371 [13 [22] ■ In [H] : it is proved that the reduced equations of motion 
of the homogeneous ball are already Hamiltonian. However, the general problem 
of integrability and Hamiltonization is still unsolved. 

1.1. Natural Nonholonomic Systems. Let Q be a n-dimensional Riemannian 
manifold with a nondegenerate metric •), y : Q — > M be a smooth function and 
let I? be a nonintegrable (n — fc)-dimensional distribution of the tangent bundle TQ. 
A smooth path q{t) G Q, t G A is called admissible (or allowed by constraints) if 
the velocity q{t) belongs to 'Dq(t) for all t E A. Let q = {qi, . . . ,qn) be some local 
coordinates on Q in which the constraints arc written in the form 
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where are independent l-forms. The admissible path q{t) is a motion of the 
natural mechanical nonholonomic system ((3,k, (or a nonholonomic geodesic 

for = 0) if it satisfies the Lagrange-d'Alembert equations 

Here the Lagrange multipliers Aj are chosen such that the solutions q{t) satisfy 
constraints (HI) and the Lagrangian is given by the difference of the kinetic and 
potential energy: L{q,q) = ^J2ij '^^j'ii'ij ~ ^il)- The expression Y,j=i^j^Hl)i 
represents the reaction forces of the constraints ([T]). 

Applying the Legendre transformation pi = dL/dqi = i^ijQj one can also 
write the Lagrange-d'Alembert equations as a first-order system on the submanifold 
M = of the cotangent bundle T*Q: 

, , . OH . OH ^--\ , , / X 

(3) " = "a^ + Z^^-j" * = l:---7": 

where the Hamiltonian is H{q,p) = ^^ij K^^'PiPj + Viq)- As for Hamiltonian 
systems, it is a first integral of the system. 

1.2. Symmetries, Chaplygin Reduction and Hamiltonization. Suppose that 
a Lie group K acts by isometrics on (Q, n) preserving the potential function V (the 
Lagrangian L is K- invariant) and let be the vector field on Q associated to the 
action of one-parameter subgroup exp(t^), ^ € t = Lie{K). The following version 
of the Noether theorem holds (see [D [3]): if is a section of the distribution V 
then 

In other words, if : T*Q — >■ t* is the momentum mapping of the X-action with 
respect to the canonical symplectic structure on T*Q, then ^k{0 conserved 
along the flow of 

On the other side, suppose that Q has a principal bundle structure n : Q ^ Q / K 
and that P is a X-invariant collection of horizontal spaces of a principal connection, 

(5) T,Q = v,(S)iq, «5 = UQ(g)|^ee}, g e Q. 

Then (Q,^, V,2?) is called a K- Chaplygin system. The system ^ is ii'-invariant 
and reduces to the tangent bundle T{Q/K) = V/K with the reduced Lagrangian 
Lred induced from L\xi- 

Let Hred be a natural mechanical Hamiltonian, the Legendre transformation of 
Lred- The reduced vector field Xred on the cotangent bundle T*(Q/K) can be 
written in the almost Hamiltonian form 

(6) ix^A^ + ^)=dHred. 
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where i7 is the canonical symplectic form on T*{Q/K), ^ is a semi-basic form 
depending of the momentum mapping $x and the curvature of the connection 2? 
(for the details see [H H H [29]). 

In some cases the equations ([2]), i.e, ([3]) have a rather strong property - an 
invariant measure (e.g, see [U |3I])- Within the class of iiT-Chaplygin systems, 
the existence of an invariant measure is closely related with their reduction to a 
Hamiltonian form. 

Suppose that the form f] + S is conformally symplectic d{J\f{^ + S)) = (it is 
assumed that 7\A is a function on Q/K). In this case the system © has an invariant 
measure Af'^~^il.'^, d = dim{Q/K) and after a time rescaling dr = J\fdt it becomes 
the Hamiltonian system with respect to the form Af{il + S) . For d — 2 the above 
statement can be inverted: an existence of an invariant measure implies that the 
nonholonomic form + S is conformally symplectic, see [9l [28l [TBI [3 [Ml [11]. The 
conformal factor J\f is called the Chaplygin reducing multiplier. 

Nonholonomic systems on unimodular Lie groups with right-invariant constraints 
and left-invariant metrics, so called LR systems, always have an invariant measure 
|30) . A nontrivial example of a nonholonomic LR system on the group SO(n) 
in- dimensional Veselova problem), which can be regarded also as a SO{n — 1)- 
Chaplygin system such that the reduced system on S"~^ = SO{n) / SO{n — 1) is 
Hamiltonian after a time rescaling, is given in |16[ (see also Section 5). 

The Chaplygin-type reduction and a (partial) Hamiltonization can be performed 
also for a class of _fr- invariant noholonomic systems {Q, n, V, 2?), where the condition 
([SJ is not satisfied on some iiT-invariant subvariety S C Q (see [T7]). 

1.3. Chaplygin Sphere and Reduction of Internal Symmetries. The n- 

dimensional Chaplygin sphere describes the rolling without slipping of an n-di- 
mensional balanced ball on an (n — l)-dimensional hyperspace TL in R" ([Il]i see 
Section 2 below). This is an R"~-'^-Chaplygin system: the kinetic energy and the 
nonholonomic distribution T) are invariant with respect to the translations of the 
ball over the hyperplane 7i. After R"~ ^-reduction it becomes the almost Hamil- 
tonian system ^ on the cotangent bundle of the orthogonal group SO{n), 

R"-i — >V ^ T{SO{n) X R"-i) 

(7) I 

P/R"-i ^ TSO{n) = T*SO(n) . 

The system is additionally invariant with respect to the SO{n — l)-action - 
rotations of the ball around the vertical vector F. The associated vector fields 
Cso(n)xR"-i ai'e sections of the connection ([7|) and we have Noether integrals ^ 
that descend to the conservation law $ = of the reduced flow. Here 

(8) $ : T*SO{n) so{n - 1)* 

is the equivariant momentum mapping of the SO{n — l)-action with respect to the 
canonical form 51 on T*SO{n). 
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However, $ is not the momentum mapping with respect to the nonholonomic 
form f7 + S. Recently, Hochgerner and Garcia-Naranjo proved that the form S 
can be truncated to the form S, such that $ is the momentum mapping of the 
5'0(n — l)-action on (T* SO{n),n + E) [22]- Moreover, the reduced system is almost 
Hamiltonian with respect to + 2 as well: ixr^di^ + S) = dHred- 

As a result, following the lines of the usual symplectic reduction, we can use the 
momentum mapping $ to reduce the system to the almost Hamiltonian system on 
(Af^,w^), where M,, = ^^^{ri)/SO{n — 1)^, SO{n — 1),, is the coadjoint isotropy 
group of 77 G so(n — 1)*: 

(9) ^x^^^^'^v = dH'^ed 

(see dH). Now H^^^ is the induced Hamiltonian function on M,,. So, the Chaplygin 
multiplier method is still applicable. In particular, if the ball is homogeneous, the 
reduced forms are closed and the reduced systems ([9]) are Hamiltonian without 
a time reparametrization. 

Let On be the coadjoint orbit of 77. The reduced space is a O^-bundle 
over T*S'"-i = T*{SO{n)/SO{n - 1)) that can be seen as a submanifold of the 
SO{n — l)-reduced space so{n)* x S"'^^: 

On Mn = <^'\On)/SO{n - 1) C {T* SO{n)) / SO{n - 1) 

(10) ; II 

T*S so{n)* X 5"'-! 

We shall consider the simplest but still very interesting and nontrivial case, 
when 77 = 0. Then the manifold AIq is diffeomorphic to the cotangent bundle of the 
sphere S"^^ and the reduced form wq is a semi- basic perturbation of the canonical 
symplectic form oj of T*S"'^^. 

For the sake of simplicity, denote wq, iJ^ed' -^red' ^"7 w, iJ, X, respectively. 

1.4. Outline and Results of the Paper. In Section 2, we recall the equations 
of motion of the Chaplygin sphere. The reduction of the system to the cotangent 
bundle of the sphere T*S'"'^^, for a zero value of the 50(n— l)-momentum mapping 
$ is described in Section 3. 

The calculation of an invariant measure as well as the time reparametrization 
dr ~ J\fdt and the reduction of the system to the Hamiltonian form for a specific 
choice of an inertia operator / of the ball is given in Section 4. On the level of 
forms, this means that the form w is conformally symplectic: d(A/'w) = 0. The 
description of the Hamiltonization is given in redundant variables, by the use of a 
Dirac bracket. 

We show that the obtained Hamiltonian system is an integrable geodesic flow. 
Moreover, as in the 3-dimensional case [12], the reduced system is closely related to 
the associated nonholonomic Veselova problem (see Section 5). Namely, the reduced 
Veselova problem and the reduced Chaplygin sphere problem share the same toric 
foliation of T*S'"'-\ 
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In the S-dimensional case, the group 5*0(2) is Abehan and ah reduced spaces 
Mj^ are diffeomorphic to T*S^. After a remarkable change of variables, Chaplygin 
transformed the problem to the case 77 = [8 . Since for n > 3 and rj Q the coad- 
joint orbits Ori are nontrivial, some additional efforts are needed for understanding 
the complete dynamics of the ball and it rest still unsolved. 

2. Chaplygin Sphere 

2.1. Kinematics. Following [TU [TS], consider the Chaplygin-sphere problem of 
rolling without slipping of an n-dimensional balanced ball (the mass center C coin- 
cides with the geometrical center) of radius p on an (n — l)-dimensional hyperspace 
H in K". For the configuration space we take the direct product of Lie groups 
SO{n) and M", where g G SO{n) is the rotation matrix of the sphere (mapping a 
frame attached to the body to the space frame) and r G R" is the position vector 
of its center C (in the space frame). For a trajectory {g{t),r{t)) define angular 
velocities of the sphere in the moving and the fixed frame, and the velocity in the 
fixed frame by 

In what follows we identify so{n) ^ so{n)* by an invariant scalar product 

(11) {X,Y) = -^triXY). 

Let / : so(n) so{n)* = so{n) be the inertia tensor and m mass of the ball. 
The Lagrangian of the system is then given by 

(12) L - iK((^, V), (cI;, V)) = + im(V, V), 

where (•, •) is the Euclidean scalar product in R". 

Let r G K" be a vertical unit vector (considered in the fixed frame) orthogonal 
to the hyperplane H and directed from TL to the center C. The condition for the 
sphere to role without slipping leads that the velocity of the contact point is equal 
to zero: 

(13) V - pflT = . 
The distribution 

V^{{g,r,Lo,V)\V = pAd,{u;)r} 

is right {SO{n) x IR")-invariant, so the Chaplygin sphere is an example of a coupled 
nonholonomic LR system on the direct product SO{n) x R" (see [25]). 

If we take the fixed orthonormal base i?i , . . . , En such that F = i?„ , then the 
constraint (jl3p takes the form = = p^m, i = ■ ■ ,n ~ 1, fn = V„ = 0, 
where ^ {Vl,E, A E^) {X AY = X ®Y -Y ® X = XY^ - YX^, X,Y e R"). 
The last constraint is holonomic, and for the physical motion we take r„ = p. 

From now on we take SO{n) x R"^^ for the configuration space of the rolling 
sphere, where R"~^ is identified with the affine hyperplane pT + T-L. Then the 
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Chaplygin sphere is an M" ^-Chaplygin system ([7]), where the reduced Lagrangian 
reads 

1 7TI fP" 1 

(14) Lred{uJ,g) = -{Iuj.,i^) + — ^(Adg(w)r, Adg(a;)r) =: -(^^6^(5) w, w). 

Remark 1. We can also consider the rubber Chaplygin sphere, defined as a system 
([T^ . subjected to the additional right-invariant constraints $7^ = 0, 1 < i < 
j < n — 1 describing the no-twist condition at the contact point [HI [25] . 

2.2. Dynamics. From the constraints (IT3|) we find the form of reaction forces in 
the right-triviahzation in which the equations ^ become 



(15) M = -pk A r, 

(16) mV = A, 

(17) 9 = ^-9, 

(18) f = V. 



where M = Adg(/w) G so{n)* = so{n) is the ball angular momentum in the space 
and A G R" is the Lagrange multiplier. 

Differentiating the constraints ([T3|) and using (|16l) we get A = mpCtV. On the 
other hand 

(19) A A r = m/9(i7r) A r = TOP (^f2 r (g) r + r r f2^ ==mppr^(f2), 

where f) C so{n) is the hnear subspace t) = M" A F and pr|j : so{n) f), pr(,(^) = 
(^r) A r = (g) r -I- r (g) is the orthogonal projection with respect to the scalar 
product pl|) . 

Whence, ([T5|) . ([T7|) is a closed system on TSO{n), representing the Chaplygin 
reduction of the R"^ ^-symmetry. Now we need to write it in the left trivialization 
of TSO{n). 

Let 7 = 5^^r be the vertical vector in the frame attached to the ball. Then 

(20) Adg-i([)) =M"A7=:f]'^. 
From the identity 

(21) w = Adg-i(0) 

and the relations ([T9|) and pr,,^ (0 = (C'7)^7 = '?7®7 + 7®7Cwe get 
Ilu = [lu!, uj] — TOp^(cj 7(g)7-|-7(g)7cj). 
Let us denote mp^ by D and let 

(22) ]<. = lu + D prfi^ u = luj + D{uj 7(g)7-|-7(g)7w) G so{n)* = so{n) 

be the angular momentum of the ball relative to the contact point (see [2]). Note 
that k = Kred{9)i^, where the reduced metric Kred{9) is defined by ^T^ . 
By using the Poisson equation 

(23) 7 = — CL17 
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it easily follows ^ {oj^ 07 + 7® 70;) = W7 (g) 7 + 7 (g) "fiu + [uj"f (g) 7 + 7 (g) "fuj, uj] . 

Therefore, the reduced Chaplygin sphere equations, in variables (k, g) of the 
cotangent bundle T*SO{n) (or in variables of the tangent bundle TSO{n)) 

are given by 

(24) k=[k,L.], 

(25) g^g-uj, 

while the reduced kinetic energy is HredO^,g) — ^{''^redid) ^j^) = hC^T^O^))- 

Let J7 be the canonical symplectic structure on T*SO{n), d — dim 50(71). It 
follows from |14[ [T5] that the reduced flow on T*SO{n) has an invariant measure 

(26) gh)\^=g-.r = 1/ Vdeti^redig)) ^'^ = 1/ ^det{I + Dwi,.)U=g-ir 

The system is additionally left 5'0(n — l)-invariant where the action of S0{n — 1) 
is given by the rotations around the vertical vector F. The closed system ([23|). (|24|) 
in coordinates (k, 7) represents the reduction of SO{n — l)-symmetry to 

(27) so{n)* X S*"-! ^ {T*SO{n))/SO{n ~ 1). 
The volume form (1261) descends to the invariant measure 



(28) f?(7) Q.oinr ^ •^s-i = l/y'det(/ + Z?pr^-,)l?,<,(„). A f^s^-i , 

where fiso{ny and fig"-'^ are standard volume forms on so(n)*(k) and S'"~^(7), 
respectively (see [TH [T5] ). 

2.3. Classical Chaplygin Sphere. In the case n = 3, under the isomorphism 
between R'^ and so(3) 

-X3 

(29) X = (Xi,X2,X3)^-~>X= I 




-X2 Xi 

from ((M)) and we obtain the classical Chaplygin's ball equations 

(30) ^k^kxw, ^7 = 7xw, 

where k = Id + Dlj — D{uj,^)^ and / is the inertia operator of the ball. In the 
space (k, 7) the density of an invariant measure ()28p is equal to 

(31) q(^) ^ l/v/det(/ + Dl) (1 - I?(7, (/ + Dl)-^^)), 

the expression given by Chaplygin in [8 . Since the system pOI) has four integrals 

(32) Fi = (k,7), ^^2 = (7,7) = 1, F3 = i(k,tJ), ^^4 = (k,k), 

it is integrable by the Euler-Jacobi theorem: the phase space is almost every- 
where foliated by invariant tori with quasi-periodic, non-uniform motion [T]. The 
integration in [5] is divided into the two steps. Firstly, equations pop are solved in 
the case the area integral Fi is zero, using elliptic coordinates on the Poisson sphere 
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F2 = 1. Then, after an ingenious linear change of variables (k,7) 1 — > (ki,7i), the 
problem transforms to the zero area case. 

3. Reduced System in Redundant Coordinates 
3.1. Reduction to T*S''-'^. From we have 



^iP'^soin-i)^) = ^(pr^o(n-i) Adg(/a;)) 

(33) = ^(pr^„(„_i) Adg(/a; + Dpr(,-,a;)) 

where so{n — 1) C so{n) is orthogonal complement to t) = M" A F with respect to 
the scalar product (fTT|) . 

The integral (l33l) is actually the momentum mapping ([8|) of the left SO{n — 1)- 
action. For ri = 3 we have the classical area integral Fi — (k,7) = (/a;, 7). 

So we can pass to the reduced system Q on = ^~^{ri)/SO{n — 1)^ ([H], see 
Introduction). We shall consider the simplest but still very interesting case, when 
we fix the value of the momentum mapping $ to be zero 

(34) prso(„_i) Adg(/w) = pr,„(„_i) Ad^ k = pr,o(„-i)7 (I^^) = Wso(n-i)-< k = 0. 

Here so{n - 1)"^ := Ad^-i so{n - 1) = (E" A 7)^ = (f)'')^. 

Whence, both k and lu belong to the subspace ((20)) . Now, let us introduce new 
variables p, £, E M" orthogonal to 7 

(35) (7,p) = (7,^) = 0, 
such that 

(36) k = 7Ap, a; = 7-^(7 AC). 
Lemma 1. The variables p and ^ are related via 

(37) P = ^-DI-'{jA^)j 

Proof. The proof directly follows from the definition k = /cj + D{{ujj) A 7) and 
relations (l36l). □ 



From ([57)) . under the conditions ([55)) . the variable ^ can be uniquely expressed 
via p and 7. 

Note that the coordinates (7,p) can be considered as redundant coordinates of 
the cotangent bundle of the sphere T*S"'~^ realized as a subvariety of K.^" defined 
by constraints 



(38) 



</'l = (7,7) = 1, (1)2 = {1,P) = 0- 
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Theorem 2. The reduced Chaply gin- sphere problem onT*S" ^ = $ ^{0)/ SO{n — 
1) is described by the equations 

(39) 7 = X^(7,P) = -^7 = A e(7,P))7 

(40) p = Xp{-f,p) = -ujp = -r\"f A C(7,P))P 

Proof. The mapping (7,p) >■ (k = 7 Ap, 7) realizes 7^*5"~i as a submanifold 
of ((27)) (see diagram ([TO|V The equation ((39|) foUows directly from the Poisson 
equation On the other hand, from the equation (|24p we get 

7Ap + 7A|j = [7AP, w] 

—OJ^p^ — p{—u)^Y' + 7 A p = 7p'^a; — p^'^LO — (jo^p^ + Lop^'^ 
'y Ap = up'y + 'yp u! ~ [up) A 7 
=^ p = — + A7. 
The multiplier A is equal to zero. Indeed, from (1381) we have 



= (7, P) + (7, P) = (-^^7, P) + (7, -'^P) + A(7, 7) = A = 0. 



□ 



Note that the reduced Hamiltonian 

(41) H{^,p) = i(k,c.) = i(7Ap,/-i(7 Ae(7,P))) 

(which is now unique only on the subvariety (1381) ) as well as the system (p9|. (|40)) 
itself, is defined on 

(42) = M^n \ {7 = 0}. 

Also considered on R^", the extended system ((39|) . (l40t preserves the functions 
01, 02, the Hamiltonian (j4ip and the reduced momentum 

(43) K{-f,p) = (7 Ap,7 Ap) = {j,l){P,P) - (7,P)^- 

3.2. Chaplygin Reducing Multiplier. At the points of the vector field 

X = {Xy,Xp) of the system (15^ . (HO]) can be written in the almost Hamiltonian 
form ix(w) = dH, where the form w is a non-degenerate 2-form on T*S"^^, a 
semi-basic perturbation of the canonical symplectic form 

(44) uj ^ dpi A d-fi -\ + dpn A d-fn |t*s"-i 

(see [12]) ■ 

Let w be an almost symplectic form, i.e., a nondegenerate 2-form on an even 
dimensional manifold M. For an almost Hamiltonian flow i = X, ixw = dH, the 
Chaplygin multiplier is a nonvanishing function Af such that cu = TVw is closed. 
Since i-^ui — dH, X — jfX, applying the time substitution dr = JVdt, the system 
X = X becomes the Hamiltonian system -^x = X with respect to the symplectic 
form u} [551 13 Uni m] • More generally, J\f is the Chaplygin multiplier if there exist a 
2-form w such that ijfW = and w = A/'(w — w) is symplectic (see [11]). Then, as 
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above, the system i — X becomes the Hamiltonian system -^x = X with respect 
to the symplectic form uj. 

Alternatively, a transparent and classical way to introduce the Chaplygin reduc- 
ing multiplier for our system is as follows (e.g., see Section 3 in [H]). Let Af{'y) be a 
difFerentiable nonvanishing positive function in a neighborhood of S'^~^. Consider 
the coordinate transformation 

(7,p) I — ^ (7,_p), p^Mp 

defined in some neighborhood of and the new symplectic form 

(45) w = dpi A H h dpn A d'y„\T*s^-^ 

= Afuj + pi dAf A dji ^ hPndA/" A d7" Ir^sr.-!. 

Then TV is a Chaplygin multiplier for the reduced system if the equations (15^ , 
(Uni) in the new time dr = N{q)dt becomes Hamiltonian with respect to the form 
Co. If A/" is a Chaplygin multiplier then from the Liouville theorem we have 

(46) /:^(c;;"-i) = o 4=> /:x(A/'""^a;"-i) = 0, 

i.e., the original system has the invariant measure with density A/'(7)"^^. Further, 
the form w reads 



(47) w = w +pi dlnA/" A d7i H V pn din M h d-y" 



T'S" 



3.3. Homogeneous Sphere. It is proved in [22] that the reduced equations of 
motion ([9]) of the homogeneous ball are already Hamiltonian, for any value of the 
50(71 — l)-momentum mapping. This interesting result, for $ = 77 = can be easily 
derived from Theorem 2. 

Suppose the inertia operator / equals si (multiplication by a constant s > 0). 
Then the equation ([57]) . under the conditions ([55]) . gives ^ = sp/{s + D). The 
reduced system (15^ . PO]) takes the form 

(48) ^-475^' 

s + U s + D 

representing the geodesic flow of the standard S'0(7T,)-invariant metric of the sphere 
multiplied by s + D. Note that in this case the angular velocity 

w = -(7 AC) = — ^:(7Ap) 
s s + D 

is constant along the flow of (US]). Actually, the angular velocity w is constant for 
the rolling of the homogeneous ball for any value of SO{n — l)-momentum mapping. 
Namely, substituting M = sil into the equations (|T5]) and (fT9)) we obtain 

pTt^{sn + Dn) =0, pr(,_L(sf2) = 0, 

which implies i7 = 0. Further, from p^ . (ITOl) . (1^ we get di — V — (see also 
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4. HAMILTONIZATION 

In this section we shall perform the Hamiltonization of the reduced Chaplygin 
sphere ([M]) , PU]) for the inertia operator defined on the base Ei A Ej via 



(49) I{E, A E^) = "'"'"^ E, AEj, l<i<j< n, 

D — Giai 



where < aiGj < D, 1 < i, j < n. 

The form of the inertia operator as well as the form of the Chaplygin multiplier 
below is motivated by the corresponding formulas in the problem of motion of the 
n-dimensional Veselova problem as well as the rubber Chaplygin ball given in [TB] 
and [25], respectively. 

Let A = diag(ai, . . . ,a„). 

In the tree-dimensional case the operator defines a generic rigid body inertia 
tensor /. Indeed, using the isomorphism ([2^ . we get 

(50) / = diag(/i,/2,/3) 

where Ii = a2a^D/{D - a2a^), h = a,^aiE> / {D - a^ai), I3 a2a3D/{D - a2a3). 

Conversely, given a generic inertia tensor (j50|) (one can always assume that the 
axes of the frame attached to the ball are principal axes of inertia), the matrix 
A = diag(ai, a2, as) is determined via 



(51) a, = ^hhhDih + D)/h^{h + D){l2 + D){h + D), i = 1,2,3. 

Remark 2. In general, for n > 4, the operator (1491) is not a physical inertia operator 
of a multidimensional rigid body (see [U). However, by taking conditions 

(52) ai = a2 = ■ ■ ■ — a„-i ^ a„. 

and 2anD > aia„ + aiD, we get the operator luj = Juj + ujj, where 

J _ r J T T\ 7 _ "1^ 7 _ aia„D ajD 

2[U ~ af) D ~ aittn 2(D — a\) 

representing a SO(n — l)-symmetric rigid body {multidimensional Lagrange case 
[2) with a mass tensor J. 



Theorem 3. The extended reduced Chaplygin sphere equations (j39[) . (|40p . defined 
by the inertia tensor (j49[) . read 



.t-o^ • 1 (P'7) <-i , h,Ap) (7,7) 

^=dP- D{^,A~W) DH-y,A'^-yf - DH^,A~^^) ^P^ 

(^A\ (P,^~W) (P,-P) .-1 {P, Ap) jp, 7) 

^ ^ £1(7,^-17)^ 15(7,^-17) £12(7,^-17)^ D2(^^^-l^) 

Proof. From the definition (1491) . the angular velocity is given by 

(55) uj = r\j A£.) = A-^j AA-^^^ A£_. 
Now, the equation (|37l) . under the conditions (|35p . can be solved 

(56) ^= Di^X'^) (^P-(P^^^h)- 
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Thus Lo = ^p- 7 A Ap/D) /L'(7, ^^^7) and simply follows from 

(EH), dini). □ 

4.1. Invariant Measure. The canonical volume form Q on M^" induces the vol- 
ume form (7 on 1^ ]]j2ii (e.g., see paragraph 3.6, Ch. 1 [T]). A simple 
description of fi, in terms of the restricted symplectic structure (j44|) is as follows. 

Consider the standard spherical coordinates {9,r) = (6*1, . . . ,f?„_i,r) on R"(7) 
and the corresponding canonical momenta (tts, tt^) = (tti, . . . , 7r„_i, tt^) on R^"(7,p) 
with respect to the canonical symplectic form: 

dpi A d'^i + ■ • • + dpn A d7n — dni A (i6'i + ■ • ■ + c?7r„_i A dOn-i + rfTTr A dr. 

Then the volume form fi can be represented as 

(57) i7 rfpi A ^71 A • • • A dpn A d-yn = (^tti A d^i A • • • A d7r„_i A dOn-i) A dp,. A dr, 



where r = VvT/t) ^i^d = (7,p)/Vw/7)- The coordinates {9,ng) are canonical 
coordinates (the symplectic form (H^ equals cZtti A d^i + • • • + d7r„_i A dOn-i) and 

(T = (iTTi A A • ■ • A dvTn-i A d9n-l 

is the canonical volume form on the cotangent bundle T*S^~^, naturally extended 
to R2». 

Proposition 4. The reduced Chaplygin system (|53|). (|54p on T*S"~^ possesses an 
invariant measure 

(58) p{j)a ^{A-'j,^r^^~'y^a. 
Proof. The divergence of the vector field X in R^" is 

(59) div(X) =±(P^ + ^)=in-2)( ^^' -'"f ) ^ + 
where 

^ = f 2(7,7) tjA-^ tr^ \ 

^£'(7,^-17)2 ^ D2(^^ ^-1^)2 i:)(-^,yl-l7) £12(7,^-17)7 



Whence, on the invariant submanifold (/)2 — vr^ = 0, in view of (|53p. we get 



i=l 



97i (7,^~V) 



In other words, the density p{'j) satisfies the Liouville equation 

„„, .M.x,.|:,.|i..|:(||.|).o 

on the manifold (j)2 ^ T^r ~ 0. 

On the other side, from (|57l) we obtain 



(61) Lx{fJ-f^) = ^xifJ'Cr) A (iTTr A dr + fia A CxidiTr A dr). 
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Since the functions (jji, 02 are invariants of the vector field X, the Lie derivatives 
CxdiTr and Cxdr equal zero. Further, (pO)) implies that the left hand side of (pT|) 
is also equal to zero on the invariant sub variety 02 = i^r =0. Thus we conclude 

Cx{^^o)\T'S^-^ = 

as required. □ 

Remark 3. The reduced vector field ([9]) has an invariant measure for $ = r; 7^ as 
well. Namely, the SO{n — l)-reduced system (IM)) preserve the volume form 

(PS)) . Then the restriction of the flow to the invariant manifold (see diagram 
pop ) preserves the induced volume form (e.g., see [T]). In this sense. Proposition 
H] is equivalent to the proportionality of the densities of measures ([25]) and ([55)1 
(compare with Theorem 5.1 [16 ). In particular, for 71 = 3 the density /i(7) after 
the transformation (j5ip . up to a multiplication by a constant, takes the form (I3ip . 

4.2. Time Reparametrization. The reduced Hamiltonian (|^T|) read 

(62) g(7,P) = 2j^(^ Vi7) ^^^^'^''^^^~ 

According to the constraints instead of (p^ we can use the Hamiltonian 

function 

(63) H{^,p) = ^jj2(^|^-i^) ^"'7)(P,P) - (P, Ap)) 

As follows from Proposition 2] and the relation pS)) . if the reduced Chaplygin 
system on T*S"^^ is transformable to a Hamiltonian form by a time reparameter- 
ization, then the corresponding reducing multiplier J\f should be proportional to 
1/V(7,^"^7). 

Theorem 5. Under the time substitution 

(64) dr =Mdt = , ^ dt 

and an appropriate change of momenta 

(65) (7,P)' — > {l,p), P= ^ n-^, , , P 

L»V(7, A-I7) 

t/ie reduced system (|53p . (|54p becomes a Hamiltonian system describing a geodesic 
flow on X*S"~^ with the Hamiltonian 

(66) 77(7,p) = i {Dij,A-'j){p,p) - (p,Ap)) . 

Proof. Consider the cotangent bundle T*S"^^ realized as a submanifold of M^" 
given by 

(67) V'i = (7,7) = l, ^2 = (7,P)-0. 

The canonical Poisson bracket on 7^*5"-i with respect to the symplectic form 
(|45|) can be described by the use of the Dirac bracket (see ^T0 \ [26 l ITj): 

{F, G}d = {F, G} - ({F, ^i}{G, V2} - {F, ^2}{G, ViD/iV'i, ^^2}, 
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where 

^ V dpi dpi dji ) 

Considered on M^", the bracket {•, - jd is degenerate and has two Casimir func- 
tions ■f/'i and -02- The symplectic leaf given by (p7)) is exactly the cotangent bundle 
rp*gn~i endowed with the canonical symplectic form. 

Under the mapping (|65l) . the Hamiltonian (1551) transforms to With the 

above notation, the geodesic flow defined by Hamiltonian function (1661) , in the time 
T, is the restriction to (l67t of 



li = -i-lt ^ {li,H}d, p- = -^Pj = i = l,...,n. 

ClT dT 

It is convenient to find equations (j68| using the Lagrange multipliers (see [26l[Tj). 
Introduce 

H* = H - X^i - Hij2- 
The equations (l68l) are then given by 

, dH* dH ..^ , 

7 = ^ ~dp ^ ^' - Ap- fij, 

OH'^ OH 

p' = — — = -_ + 2A7 + = -D(p,p)v4-i7 + 2A7 + /^p 
07 07 

where the multipliers A and /i are determined from the condition that the constraint 
functions ipi and "02 are integrals of the motion. 
Straightforward calculations yield 

_ {Ap,p) _ D{-f,A-^-f){p,j) - {Ap,j) 



2(7,7)' (7,7) 



and therefore 



(69) y = D(A-S,7)p-Ap+ 7 7 



(7,Ap) L>(7,A-i7)(p,7) 



(7,7) (7,7) 



-7, 



(70) p=-D{p,p)A 7+- ^7 "-7 rP+ 7 ^ P- 

(7,7) (7,7) (7,7) 

In the time i, inverting the mapping (j65p . the equation (j69p takes the form 
7-C\/(7,^-^7) ^ ^ 



i?V(7,^-^7) 

Z?(A-7,7)p-Ap+(^'^^)-^ i?(7,A-7)(p,7), 



(7,7) (7,7) 
i.e., 

. ^ J_ 1 . (7,^P) (P,7) 
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which coincides with (|53p at the points of T*S" ^. Further, 



(72) 




dV(7,^-h) V(7,^-h); "(7,^-i7)^ 

Finally, substituting p = A/p into the right hand side of ([70]). combining with 
(ini) and (Hal), we get 

(73) p = ^^'P) + ^1^^^^ 7 

^ ^ ^ ^(7,^-^7) D^{j,A~^j){^,^)^ 

D{j,A-^-ff D^j,A~^j)^^' 
As above, the equations and ((75)) are different, but they coincide on the 
invariant manifold (pi — — 1, (f>2 — V'2 = 0. The theorem is proved. □ 

Remark 4. The link between the Dirac bracket and the Lagrange multiplier ap- 
proach can be expressed via 

{g>2} {g,V'l} 

Also, note that the reduced almost symplectic form (|47p is given by: 

Pitt'^^j 

w = 2^ A dj, - - — djj A dj, |t*S"-i- 

Remark 5. During the referee process of this paper, the paper [23| appeared, where 
the Abelian O-Chaplygin systems associated to Cartan decompositions g = t © p of 
semi-simple Lie algebras are studied. They are defined on the direct product of a Lie 
group K ^ Lie{K)) endowed with a left-invariant metric with the vector space D 
(0 = [r, t] C p) endowed with the metric induced from the Killing form. Here F G p 
is fixed. As an example, taking the Cartan decomposition so(n, 1) = so(n) © K" 
of the Lie algebra so(n, 1) one gets the Chaplygin sphere problem (compare with 
the equations (73) in [IS]). Besides 0-reduction to T*K, likewise the Chaplygin 
sphere problem, the system has an internal symmetry group H <Z K (isotropy 
group of F) and admits the almost symplectic reduction with respect to the H- 
action. Hochgerner derived the equations on the parameters of the kinetic energy, 
such that the (zero momentum) reduced almost symplectic form is conformally 
symplectic. The operator ([^^ represents the solution of these equations within the 
class of diagonal operators on so{n) with respect to the base Ei A Ej [23]. 

5. Integrability 



5.1. Classical Chaplygin Sphere and the Veselova Problem. The Veselova 
problem describes the motion of a rigid body about a fixed point subject to the 
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nonholonomic constraint 

(74) (w,7) = 0, 

where w is the vector of the angular velocity in the body frame and 7 is a repre- 
sentation of a unit vector fixed in a space, relative to the body frame [30] . The 
equations of motion in the moving frame have the form 

(75) — Iw = Iw X w + X-y, —7 = 7 x w, 
^ ' dt ' dt' ' ' 

where I is the inertia tensor of the rigid body and A is a Lagrange multiplier chosen 

such that w(i) satisfies the constraint ([74]) . 

f7R\ \ (Xw X w,X-^7) 
76) A = . 

[I ^7,7) 

Here we suppose that all eigenvalues of I are greater then 1. 

Equations (iTSl) . ((76)) also define a dynamical system on the whole space R^(w, 7), 
and the constraint function fi — (w, 7) appears as its first integral. The system 
has an invariant measure with density -y/ 7). Following [121, by introducing 
K — Iw — (X^w, 7)7, = I — I one can write system (|75|) . (|76|) in the form 

d ^ ^ ^ d ^ ^ ^ 

(77) — K = K X w, — 7 = 7 X w. 
dt dt 

Apart from /i = (w,7) = {K,"f), it always has the geometric integral /2 = 
(7, 7) = 1 and two other independent integrals 

(78) /3 = i(i?,w)-i(i?,7)(Z0w,7), h^{K,K). 

On the constraint subvariety (|7^. these functions reduce to the energy integral 
i(Jw,w) and (/w,/w) - (/w,7)2 (see [5D]1. 

By the Euler- Jacobi theorem [T] , the above system is solvable by quadratures on 
the whole space R^. For /i = the system was integrated by Veselova (e.g., one can 
find the motion using the isomorphism with a celebrated Neumann system [30]). 
Next, as was shown in J2], the restriction of the extended Veselova system ([75| . 
([7S|l onto the level variety /i = ci (ci 7^ 0) can be reduced to this system on the level 
/i = by a linear change of variables (if, 7) 1 — > {Ki,'yi) and an appropriate time 
reparametization. This linear change was found by using a relation of the Veselova 
system with the Chaplygin sphere problem, which we are going to describe now. 

Define the operator / and vector lj by: 

(79) I=I + Dr^, w = -luj I = D{I-iy\ w = --^(Z-I)w. 
Now we can state the following remarkable correspondence: 

Theorem 6. (Fedorov |T^) The invariant tori fi — ci, f2 = I, fs = c^, ~ C4 of 
the Veselova problem ([75t . ([75)) . via ([7^ transform to the invariant tori Fi = ci, 
F2 — 1, F3 — C3, F4 — C4 of the Chaplygin sphere problem ([50)) . 
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Let US mention that there are two interesting isomorphisms between the Chap- 
lygin sphere problem ([50)) with Fi = (k, 7)^0 and the Clebsh case of the Kirchoffs 
equations of a rigid body motion in an ideal fluid, with a zero area integral. The 
first one is described in [13. and the other one is given recently in [5]. 

5.2. Veselova Problem on SO{n). It appears that the analogue of Theorem [5] 
can be formulated for an arbitrary dimension n and a zero value of the momentum 
([M)) . First, for a reader's sake, we shall briefly recall some definitions and results 

of [ini. 

Consider a nonholonomic LR system on SO{n) defined by the left-invariant La- 
grangian Lx{g,g) — ^(Iw,w) = — itr(Iww) where I : so{n) — )■ so{n) is positive 
definite and the right- invariant distribution T>r on TSO{n) whose restriction to the 
algebra so{n) is given by c) = span{i?i /\Ej | i = 1, . . . , r, j = 1, . . . , n}. This implies 
the constraints 

(80) (w, Adg-i {E, A Ej)) = (w, e, A e^) = 0, n-~r + l<i<j<n. 

Here w{t) = ■g{t) G so{n) and ei = (en, . . . , ei„)'^, . . . , e„ = (e„i, . . . , e,™)"^ 
is the orthogonal frame of unit vectors fixed in the space and regarded in the moving 
frame {Ei =5-61, . . . , En = g • e„, where Ei = (1,0, .. . ,0)"^, . . . , En = 
(0, . . . , 0, 1)"^). They play the role of redundant coordinates on SO{n). 

The system is described by the kinematic Poisson equations 

(81) e-i = -vfei, i = 1, . . . , n, 

together with the Euler-Poincare equations with indefinite multipliers Xpq 

(82) ^ {Iw) = [Iw, w] + Xpg Bp A eg. 

n—r-\-l<p<q<n 

Since for n = 3, r = 2 the above system represents Veselova problem, we refer 
to {SO{n), Lxj'Dr) as a generalized Veselova system (see Fedorov and Kozlov pjj). 

The Lagrangian Li and the distribution are invariant with respect to the 
left SO{n — r)-action, where SO{n — r) is the subgroup of SO{n), rotations that 
leave Ei, . . . , Er invariant. Moreover, the distribution Vr can be seen as a principal 
connection of the bundle 

SO{n -r) — ^ SO{n) 
i 

Vn,r = SO{n)/SO{n - r). 

As a result, the system can naturally be regarded as a Chaplygin system and 
dynamics is reducible to the Stiefel variety Vn^r- The points of the Stiefel variety 
can be seen as matrices X = (ei, . . . , e^) (positions of the r-frame given by vectors 
(ei, . . . , e,-)). So, the tangent bundle TV{r, n) is the set of pairs {X, X) of n x r 
matrices subject to the constraints 

(83) X'^X = Ir, X'^X + X^X = 0. 
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The reduced Lagrangian takes the form Lred{X , X) ~ — |tr(I$,.$r) (see flBI). 
where the tangent bundle momentum mapping $r '■ TVn^r ~^ so{n) = so(n)* is 
given by 

<i>r{X,X) = XX'^ - XX'^ + ^X[X^X - 

= eiAeiH h A + ^ ^ [{ea,ep) - {ea,ep)] Ca A ep. 

l<a</3<r 

Introduce the n x r momentum matrix 

(84) V,,=dLred{X,X)/dX,,. 

Since the Lagrangian is degenerate in the redundant velocities Xis , from this relation 
one cannot express X in terms of (A", V) uniquely. On the other hand, the cotangent 
bundle T*V{r, n) can be realized as the set of pairs {X, V) satisfying the constraints 

(85) x^x = ir, x^v + r^x = 0. 

Under the conditions ((83|) . (|85)) . the relation ((84)) can be uniquely inverted, i.e., one 
gets X = X{X,V). Then we have (see Theorem 5.4 in [161): 

Theorem 7. ([16 ) The SO{n - r)- reduction of the Veselova problem ([5T]) . 
(1821) is the restriction to T*V{r,n) of the following system on the space {X,V): 

(86) X ^ ^^r{x,v)x, r ^ -^rix,p)v, 

where <S>riX,V) = <i>riX,X{X,r)). 

Remark 6. Here we use the opportunity to mention one correction to [161 : in 
the equation (5.21) the momentum mapping $* should read $* = Iw|-d,, = XV^ — 
VX^ . This equation was used only in the proof of Theorem 5.4 [T5]. The statement 
of the theorem itself remains to be correct. 

In particular, for r = 1, the Veselova problem is reducible to T*5"'~^. 
Let, as above, A — diag(ai, . . . , a„) and denote ^ ^ X = ei, p = V . Taking the 
special inertia operator defined by 

(87) I{E, A Ej) = ar^aj^ E, AEj, l<i<j< n, 

we have $i = 7A7 = jAAp/ (7, ^"^7) and the reduced system ([55)) becomes (here 
we replaced the matrix A from [Tn|[5n] by A^^) 

(88) 7 = -$i7= ^— ^j-y (-(p,A7)7 + (7,7)Ap), 

(89) p = p = {{p, Ap)j - {p, j)Ap) . 

Furthermore, as it follows from (THl dO] , under the time substitution ([M)) and the 
change of momenta (|65p the reduced system transforms to a Hamiltonian system 
describing an integrable geodesic flow on T*S"~^ with the Hamiltonian 

(90) Hh,p) = ^iAp,p). 
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Remark 7. The reduced Veselova system (j88| . (j89|) is trajectory equivalent to the 
geodesic flow on the eflipsoid i?"^^ — {x — (xi, . . . , a;„) G M" | (x. Aa;) 1}: the 
geodesic hnes x{t) of the ehipsoid, after the Gauss mapping ^{t) ~ Ax{t) l\A'y[t)\ 
and a time rescahng, become solutions of the reduced Veselova system, and vice 
versa (see [ID])- 

5.3. Integrability of the Reduced Chaplygin Sphere Problem. Let us sup- 
pose ai ^ CLj, i ^ j- As in the three-dimensional case [12] . we have 

Theorem 8. (i) The geodesic flow (|69|). ()70|) is completely integrable onT*S"~^. 

(ii) The zero SO{n — l)-momentum reduced multidimensional nonholonomic 
Chaplygin sphere problem (|53p . (j54p defined by inertia operator (1491) and the re- 
duced Veselova problem (1551) . ([5^ defined by inertia operator (|57| have the same 
invariant toric foliation of T* S"~^ . 

(iii) Let T be a regular, {n — 1) -dimensional invariant torus. Then there exist 
angle coordinates ipi, . . . , ipn-i on T in which both problems simultaneously take 
the form 



with frequencies uj\, . . . , i — 1,2, respectively. 

Proof. In what follows, we restrict our considerations to T*S^~^ . The momentum 
integral (j43| in variables (7,73) becomes 

(91) K{^,p)=D\A-^^,^){p,p) 

and the Hamiltonian (j66p can be written in the form 



H{-i,p) = ^K{j,p) - ^V.{-i,p). 



Since ^ is the integral of the geodesic flow on T*^""^ we have 

{H, K]d = {H, n}d - {n, K}d = 0. 

Consider the spheroconical coordinates (Ai, . . . , A„_i) (ai < Ai < 02 < • ■ • < 
A„_i < a„) on S"^^ defined by the relations 

(92) ^,^ {a.~X,).--{a.-X . = 1,...,. 

(see |26|). Let (/ii, . . . , /i„_i) be the canonical momenta on the cotangent bundle 
with respect to the form (P5|) 

a) — dpi A d7i + • ■ • + dpn A d7„|x. 5,1-1 = rf^i A dXi + • ■ • + dfin-i A dXn-i. 
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Then, according to [26], [6] and [16], respectively, we have: 

/- . (Afc - oi) ■ • ■ (Afe - o„) 2 
{p,p) = -4 2^ = — Mfe, 

, \ A1A2 . . . A„_i 

(7, ^ 7) = , 

aia2 . . . a„ 

^ n.^fc(A.-A.) 

Therefore, the Hamiltonian (P(7| has the Stackel form in spheroconical variables 
and the geodesic flow on determined by H is completely integrable (see 

[16|). We have Poisson commuting, quadratic in momenta integrals Fi, . . . , F^-i 
(e.g., see [T]). One can prove that functions Fi commute with iJ using the direct 
calculations in elliptic coordinates. 

Alternatively, note the geodesic flow of H, over a generic invariant torus T (level 
set of Fi, . . . ,-F'„_i) is quasi-periodic with non-resonant frequencies (for example 
this follows from Remark[7|). Thus, since {'H, K^^ = and the integral trajectories 
are dense on T, K is also constant along the Hamiltonian flows of Fi over T. Since 
we deal with analytic functions, we get that K is in involution with ^Fi, . . . , Fn-i 
on the whole T*S"~^ {K is the analogue of the classical Joachimsthal's integral of 
the geodesic flow on the ellipsoid ^M)- Further, the Hamiltonian H, as a 

linear combination of K and Poisson commutes with Fi as well. Whence, the 
system dBSj), GO]) is completely integrable on T*S"'^^. 

The last assertion of the Theorem follows from the Liouville- Arnold theorem [1] 
and the fact that the systems transform to a Hamiltonian form after the same time 
reparametrization □ 

The system is integrable even if not all a,; are distinct. For any pair of equal 
parameters ai — aj, the geodesic flow (|69| . ([70| has the additional linear integral 
fij — liPj ~ljPi - For example, let n = 4 and ai = 02 ^ 03 04. Then the complete 
set of commuting integrals is /12, /34 and H. If we have at least three equal 
parameters, the system is integrable according to the non-commutative version of 
the Liouville theorem. 



Remark 8. Note that the operators (|49| and ((87)) are related via 

DI^l + Dr^ I^D{DI-I)-\ 

In order to reobtain Fedorov's correspondence (17^ for n = 3 and fi — Fi — 0, 
instead of ([57]) one should consider the inertia operator multiplied by D. 



5.4. Lagrange Case. Consider the Lagrange case (1521) . Due to the additional 
SO{n — l)-symmetry, the geodesic flow (|69|. ([70]) has the integrals fij, 1 < i < j < 
n — I. Thus, in the original coordinates we get integrals 

(93) = (7, A-^j){j,p, - 7,pO', l<i<j<n-l. 
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In this case we do not need Hamiltonization to integrate the reduced system, 
it is aheady integrable according to the Euler-Jacobi theorem. Since the generic 
invariant manifolds given by H and integrals (|93p are two-dimensional and the 
system has an invariant measure we have [1]: 

Theorem 9. The Lagrange case of the reduced Chaplygin system ()53|) . ()54|) is 
solvable by quadratures; compact regular invariant manifolds given by functions 
(1931) and (j66p are two-dimensional tori. 

Remark 9. Although the Lagrangian (jl2p is additionally invariant with respect to 
the right SO{n — l)-action, the integrals ([M)) are not Noether's integrals. The 
reason is that the associated vector fields do not satisfy constraints (IT51) . For n = 3 
and Ii = I2, the corresponding integral of the system (|30p is = fc| — -D(7, (/ + 
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